RA signaling pathway combined with Wnt signaling pathway regulates human-induced pluripotent stem cells (hiPSCs) differentiation to sinus node-like cells

Background The source of SAN is debated among researchers. Many studies have shown that RA and Wnt signaling are involved in heart development. In this study, we investigated the role of retinoic acid (RA) and Wnt signaling in the induction of sinus node-like cells. Methods The experimental samples were divided into four groups: control group (CHIR = 0), CHIR = 3, RA + CHIR = 0 andRA + CHIR = 3. After 20 days of differentiation, Western blot, RT-qPCR, immunofluorescence and flow cytometry were performed to identify sinus node-like cells. Finally, whole-cell patch clamp technique was used to record pacing funny current and action potential (AP) in four groups. Results The best intervention method used in our experiment was RA = 0.25 µmol/L D5-D9 + CHIR = 3 µmol/L D5-D7. Results showed that CHIR can increase the expression of ISL-1 and TBX3, while RA mainly elevated Shox2. Immunofluorescence assay and flow cytometry further illustrated that combining RA with CHIR can induce sinus node-like cells (CTNT+Shox2+Nkx2.5−). Moreover, CHIR might reduce the frequency of cell beats, but in conjunction with RA could partly compensate for this side effect. Whole cell patch clamps were able to record funny current and the typical sinus node AP in the experimental group, which did not appear in the control group. Conclusions Combining RA with Wnt signaling within a specific period can induce sinus node-like cells. Supplementary Information The online version contains supplementary material available at 10.1186/s13287-022-03006-8.

Selecting RA and Wnt signaling pathways as experimental subjects is theory-based. Firstly, RA is expressed in the venous pole of the posterior second heart field and regulates the left-right asymmetry of the stages of cardiac development [9][10][11]. Secondly, a study by Protze SI showed that RA can enhance the phenotype of Open Access  13:324 biological pacemakers. However, due to the complexity of the protocol of that study, it is necessary to find a simpler method of pacemaker generation [8]. Thirdly, Ren' study illustrated that canonical Wnt5b signaling directed outlying Nkx2.5 + mesoderm into pacemaker cardiomyocytes [7]. Finally, CHIR synergizes with RA during cardiac development, while no study has elaborated their effects on the formation of a biological pacemaker. The main aim of our experiment was to determine whether combining RA and Wnt signaling pathways can help to induce sinus node-like cells and the biological mechanism underlying the induction.

Cell culture and identification Preparation before cell culture
In every six-well culture plate, 1 ml of coated medium (Aliquot Matrigel (Cat. no. 354277; Corning) + Knockout DMEM + 1% penicillin/streptomycin) was added. The liquid was quickly released at the center of the bottom of the culture plate, which was shaken immediately after adding the medium to each well, allowing the liquid to evenly spread throughout the culture plate. The whole process was performed quickly to avoid rewarming of the Matrix gel. A plate coated with Matrigel can be left in a refrigerator at 4℃ for no more than one week.

Cell passage and culture
When the cells reached 80-90% confluence, they were passaged. Versene (Cat. no. 15040066; Gibco) digestion solution was added to digest at 37℃ for 6-8 min. When the cells turned spherical and the intercellular gaps became larger, the digestion solution was aspirated, and 1 ml of subculture medium was added to each well to stop digestion. The passage ratio was generally 1:5 to 1:10. The cells in the culture plates were shaken before and after for uniform distribution. The medium was removed and replaced daily by a fresh mTeSR1 medium (Cat. no. 85850; STEMCELL Technologies). On the first day, ROCK inhibitor Y27632 (Cat. no. 72302; STEM-CELL Technologies) was added to the culture medium to facilitate cell survival.

Cell differentiation
Based on previous protocols and former experience in our laboratory, hiPSCs or hESCs were differentiated into cardiomyocytes by manipulating Wnt/β-catenin signaling [12,13]. Briefly, at the beginning of cardiac differentiation (Day 0), hiPSCs or hESCs were cultured in the basic differentiation media (containing Roswell Park Memorial Institute [RPMI]-1640 medium, B27 minus insulin (Cat. no. A189560; Gibco), L-ascorbic acid 2-phosphate sesqui-magnesium salt hydrate (Cat. no. A8960; Sigma-Aldrich)) combined with Wnt/β-catenin signaling activator CHIR99021 (Cat.no.2052; STEMCELL Technologies) for 48 h. After two days, the medium was replaced by a basic differentiation medium for 24 h. Then, the medium was replaced by a cardiac mesoderm induction medium that consisted of RPMI-1640/B27 minus insulin with Wnt pathway inhibitor IWR-1 (Cat. no. 72562; STEMCELL Technologies) for another 48 h. On day 5, the medium was removed and the cells were maintained in RPMI-1640/B27 minus insulin without any additional factors for two days. After that, the cell cultures were changed to RPMI-1640/B27 supplement with insulin (Cat. no. 17504044; Gibco) for the next few days and refreshed every other day. The cell morphology and cell pulsation in each stage of differentiation were observed and photographed by a microscope, and the markers of cells at each stage of differentiation were identified. D20 was the end observation time.

Reverse transcription-quantitative polymerase chain reaction (RT-qPCR)
Total RNA was extracted from the differentiated hiPSCs or hESCs on D20 using TRIzol ® reagent (Invitrogen). Quantitative PCR was performed to evaluate the mRNA expression of SRY-box transcription factor 2(SOX2), octamer-binding protein 4 (OCT-4), mesoderm specific T-box transcription factor Brachyury(BryT), mesoderm posterior BHLH transcription factor 1(Mesp1), cardiac troponin T(cTNT), gene encoding atrial natriuretic peptide(NPPA), sodium voltage-gated channel alpha subunit 5(SCN5A), myosin, light polypeptide 2(MYL-2), human short stature homeobox 2(Shox2), insulin gene enhancer binding protein 1(ISL-1), T-box 18(TBX18), TBX3, TBX5, hyper-polarization activated cyclic nucleotide-gated potassium channel 4(HCN4), paired-like homeodomain transcription factor 2(Pitx2) and related iron channels as well as connexion. Isolated RNA (2 µg) was converted into cDNA using a First Strand cDNA Synthesis Kit (Toyobo Life Science) as follows: 42 °C for 2 min, 37 °C for 15 min, 85 °C for 5 min, and 4 °C for 10 min. The primers used for PCR amplification were synthesized by Invitrogen Biotechnology (Thermo Fisher Scientific) and are presented in Table 1. RT-qPCR was performed using the Step-One ™ Real-Time PCR system (Life Technologies, Carlsbad, CA, USA). The reactions were then conducted using the SYBR ® Premix Ex Taq TM II (Takara Bio, Japan). Protocol was as follows: predenaturation at 95 °C for 3 min, denaturation at 95 °C for 30 s, annealing at 58 °C for 10 s, and a final extension at 72 °C for 30 s, a total of 40 cycles. Semi-log amplification curves were analyzed using the 2 −ΔΔCt comparative quantification method and the expression of each gene was normalized to β-actin. PCR analyses were repeated at least three times to verify results.

Western blot analysis
The differentiated hiPSCs on D20 were plated on 12-well culture dishes. The cells were harvested using RIPA lysis buffer (Beyotime Institute of Biotechnology, Haimen, China). Equal amounts of proteins were loaded onto a gel for sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and the separated proteins were transferred onto a nitrocellulose membrane, and incubated with primary antibodies against HCN4 (ab32675; Abcam, Cambridge, MA, USA), TBX18(ab115262; Abcam, Cambridge, MA, USA), ISL-1(abs132916; Abcam, Cambridge, MA, USA) TBX3(ab154828; Abcam, Cambridge, MA, USA), TBX5(Gene-Tex; Cat No. GTX113849), α-tubulin (GB11200; Service-bio, Wuhan, China), GAPDH antibody (Abcam; Cat. no. ab37168) and Shox2(ab55740; Abcam, Cambridge, MA, USA) overnight at 4 °C. The primary antibodies were detected by incubating the membrane with secondary antibodies for 1 h, and then enhanced chemiluminescence detection (ECL; Beyotime Institute of Biotechnology) was performed. The level of α-tubulin or GAPDH was considered to normalize the signal intensities. The Western blot assay was repeated at least three times to verify the results.

Immunostaining studies
The differentiated hiPSCs on D20 were placed on gelatin-coated coverslips in 12-well culture dishes. The cell cultures were washed with PBS and fixed with 4% paraformaldehyde, followed by permeabilization with 0.1% Triton X-100 for 15 min and blocked using bovine serum albumin solution for 30 min at room temperature. Then the cells were incubated with the primary antibodies Table 1 Polymerase chain reaction primers used in this study

Gene
Primer (  . We randomly selected three visual fields from three different cell isolates to observe the positive rates.

Patch clamp
Whole-cell electrophysiology was performed as described. The cells were plated on gelatin-coated coverslips in 24-well culture dishes. The whole cell patchclamp technique was used to record the funny current (I f current) after 2-4 days of re-plating. We selected smooth and plump cells to record the current. The impedance of the fluid-filled electrode was 5 to 8MΩ. The experiments were performed using an Axon patch-clamp amplifier 700B ( to record I f current. The Clampex program was used to analyze the sample. The sampling frequency was 10 kHz, and the filtering rate was 5 kHz. I f currents were recorded by holding the resting membrane potential at − 40 mV, then stepping to a test voltage of − 140 mV for 2 s with a step of 10 mV in each sweep.

Flow cytometry analysis
The differentiated hiPSCs or hESCs on D20 were digested with 0.25% trypsin and 0.5 mM ethylenediaminetetraacetic acid to form a single-cell suspension, washed with PBS, and centrifuged at 300 g for 5 min. After centrifugation, the cells were fixed with 4% paraformaldehyde for 30 min at room temperature, then permeabilized and blocked by 100ul blocking solution (PBS + 5% FBS + 0.2% trelaton + 0.1% Triton X-100) overnight at 4 °C. The cells were incubated with mouse anti-cTNT or mouse anti-Shox2 (Immunoway Cat no.YT6364) and rabbit anti-Nkx2.5(E1Y8H#8792; Cell Signaling Technology, Massachusetts, USA) in the dark for 1 h at room temperature, washed with PBS twice and followed by incubation with the corresponding secondary antibodies for another 30 min. Finally, the cells were resuspended in PBS before performing a flow cytometer assay (Beckman CytoFLEX).

Statistical data analysis
GraphPad Prism7 was used for data analysis. The reported data were expressed as the means ± SEM. The statistical significance of the differences between groups was determined using the one-way ANOVA. All data were subjected to formal tests for normality. Data not exhibiting a normal distribution were evaluated by nonparametric tests. P < 0.05 was considered to indicate a statistically significant difference.

Differentiation of hiPSCs into CMs
In this study, following a previously designed protocol in our laboratory, we used Wnt signaling activator CHIR and Wnt signaling inhibitor IWR-1 to induce cardiac differentiation (Fig. 1a). The pluripotency marks (SOX2 and OCT4) showed the highest expression on D0 (Figs. 1c, d).
Expression of the mesoderm-specific T-box transcription factor Brachyury (BryT) was observed on from D1 to D3 of differentiation, which indicated that the hiPSCs had undergone mesodermal induction (Fig. 1e). Mesp1, the characteristic marker of the cardiac mesoderm stage increased from D3 to D5 (Fig. 1f ). Early cardiomyocyte makers (Nkx2.5 and CTNT) were detected as early as on D7 of differentiation (Fig. 1g, h).

The optimal RA intervention time
Previous studies have found that from D3to D5, RA (HY-14649; MedChemExpress) intervention can enhance the pacemaker phenotype [8]. Therefore, our experiment initially followed this protocol to add RA = 0.25/0.5/1.0 µmol/L at days D3-D5 to induce the differentiation of hiPSCs (Fig. 1b). We found that the atrial muscle cell marker differentiated with RA at D3-D5 elevated relative to its level in the control group (CHIR = 0 group), but the levels of the pacing markers, such as Shox2 and TBX18, almost had no change (Fig. 1i-t). Several related articles about RA and cardiac embryonic development and showed that simple supplementation with RA at the embryonic stage E8.5-E14.5 could remedy the development of the SAN in RALDH2 knockout mice [14]. Moreover, a study by Ren showed that CHIR = 3 µmol/L D5-D7 can increase the generation of hiPSC-pacemaker cardiomyocytes [7]. Based on the results of those studies, we adjusted the protocol. The RA intervention period changed to RA = 0.25 µmol/L D5-D7, D7-D9, D9-D11, D5-D11, respectively, combined with CHIR = 3 µmol/L D5-D7 or not (Additional file 1: Figure  S1a). RNA results showed that both D7-D9 and D9-D11 groups could improve the expressions of pacing-related Fig. 1 Differentiation of hiPSCS into CMs and exploration of the optimal RA intervention time. Data points present as means ± SEM (n = 1). a Schematic of the experimental protocol for the differentiations of cardiomyocybytes from hiPCs with small molecule modulators of Wnt signaling. CHIR (CHIR99021) and IWRI, agonists, and inhibitors of the Wnt pathway, respectively. b Schematic of the experimental protocol for RA = 02.4/0.5/1.0umol at days D3-D5 to include hiPSCs differentiation, c-h Samples of mRNA were collected and measured from differentiating cells between day0-9 using qRT-PCR and normalized to TBP expression. i-t The4 transcriptional level of the central atrial muscle cell marker and pacemaker markers. Data points present as means ± SEM (n = 1). *p < 0.05 compared to control group (CHIR = 0 group). RT-qPCR, quantitative reverse transcriptase-polymerase chain reaction transcription factors and the pacing proteins HCN4 (Additional file 1: Figure S1b). To further determine the optimal experimental protocol of RA intervention, we expanded the RA intervention time to D5-D9, D7-D9 and D7-D11 (Additional file 1: Figure S1c). As the results showed that the mRNA levels of Shox2 and HCN4 in D5-D9 increased significantly compared to their levels in the other groups (Additional file 1: Figure S1d). Then we used RA = 0.25 µmol/L D5-D9 + CHIR = 3 µmol/L D5-D7 as the optimal RA intervention time.

The optimal RA intervention concentration
Based on the optimal intervention times described above, we further investigated the optimal intervention concentrations of RA. We used three concentration gradients of RA = 0.25, 0.5, and 1.0 µmol/L for observation. The results showed that pacing-related transcription factors and proteins were the most obvious at RA = 0.25 µmol/L group, and there was no significant effect on the expression of CTNT (Additional file 2: Figure S2a, S2b).

RA promoted the differentiation of hiPSCs into sinus node-like cells
Combined with the above experiments, we finally determined that the best intervention was RA = 0.25 µmol/L D5-D9 + CHIR = 3 µmol/L D5-D7. Then we divided the samples into four groups: CHIR = 0, CHIR = 3, RA + CHIR = 0, and RA + CHIR = 3 (Fig. 2a). The results of Fig. 2b showed that both CHIR and RA increased the mRNA level of TBX18, while CHIR had stronger influence than RA. Experimental studies illustrated that administering RA alone can increase the mRNA level of Shox2 and the pacing protein HCN4, and reduce the level of Nkx2.5 and Pitx2. The effects of the combined treatment of RA and CHIR were more prominent, and the results were statistically significant ( Fig. 2c-f ). In addition, the results of the RT-qPCR analysis demonstrated that CHIR can increase the mRNA level of Tbx3 and ISL-1, but in RA + CHIR = 0 and RA + CHIR = 3 group, the levels of both TBX3 and ISL-1 decreased (Fig. 2g, h). Surprisingly, RA could dramatically enhance the level of CTNT (Fig. 2i). The expressions of the TBX18 and ISL-1 proteins in four groups were similar (Fig. 2j, k, l). RA substantially enhanced the expression of Shox2 and slightly increased TBX3 (Fig. 2j, m, n). The HCN4 level increased in the other three groups compared to that in the group with CHIR = 0 (Fig. 2j, o).

RA rose the frequency of beats within the cardiac cell mass
The results of the flow cytometry and immunofluorescence assays further illustrated that RA combined with CHIR significantly increased the proportion of hiPSCpacemaker cardiomyocytes (CTNT + Shox2 + Nkx2.5 − ), where CHIR reduced the expression of NKx2.5, while RA mainly increased Shox2 and slightly reduced Nkx2.5 (Fig. 3). Although CHIR can increase the proportion of pacing cells, it partially decreased the frequency of cell beats, and combined with RA can compensate for this side effect (Fig. 4a, Additional file 3 : Video S1, Additional file 4: Video S2, Additional file 5: Video S3, Additional file 6: Video S4).

Effect on the expression of ionic channels
We further examined the effects of RA on the expressions of the relevant ionic channels. RA and CHIR interventions reduced the expression of sodium channels without significantly affecting potassium channels (Fig. 4b, c). Additionally, RA, along with CHIR, significantly amplified the transcription level of the calcium ion channels, cardiac sodium-calcium exchanger (NCX) and connexin 45 (Cx45), and reduced the expression of Cx43 (Fig. 4d-h), which might be the ion mechanism by which RA increases the frequency of beating.

Recorded funny current and action potential in four groups
We investigated the electrophysiological properties of the four groups of cells using the patch-clamp technique. The results showed that the funny current could be recorded in the later three groups (CHIR = 3, RA + CHIR = 0, and RA + CHIR = 3). The highest current was recorded in the RA + CHIR = 3 group, and the current decreased with the reduction of the given voltage stimulation (Fig. 5ae). The morphology AP of the control group was similar to that of working cardiomyocytes, while the typical sinus node AP was recorded in the other three groups (Fig. 5f ). Furthermore, the parameters of AP showed that AP amplitude (APA), the maximum diastolic membrane potential (MDP) level and 0-phase depolarization rate (dV/dt) in the later three groups were smaller than that in the CHIR = 0 group, while AP duration (APD 90 ) was not significantly different (Fig. 5g).

The biological mechanism by which RA promotes pacing cell differentiation
To determine the potential biological mechanism by which RA increased the proportion of pacemaker cells, we analyzed several transcription factors and miRNAs associated with the development of SAN. We found that both RA and CHIR reduced the transcription level of miR-106b, but had no significant effect on the content of miR-17-92 and miR-1 (Fig. 6a-c). Additionally, RA significantly enlarged the levels of TBX5 and BMP (Fig. 6d-e). Dynamic changes of TBX5 from D0 to D19 showed that TBX5 increased continuously from D5 till the final observation (Fig. 6f ). The results of Western blot assays demonstrated that RA significantly enhanced TBX5 expression (Fig. 6g, h). We further investigated the changes of cardiac pacing-related factors from D0 to D19 dynamically. We found that RA substantially increased Shox2 from D9 and HCN4 after D13, while CHIR influenced TBX3 after D11 and ISL-1 after D15. Both RA and CHIR reduced Nkx2.5 and slightly increased TBX18 from D9 (Fig. 7).

Discussion
Shox2 is expressed in the sinus region during mouse embryonic stage E5.5 [15]. Shox2 mutation leads to the death of the embryo at about E11.5, due to a failure of SAN [16]. The phosphorylation of Shox2 can inhibit the expression of Nkx2.5, which prevents the transformation of sinus node cells into working cardiomyocytes [17]. Moreover, cardiac-specific knockout Shox2 gene decreases the heart rate and changes gene expression patterns, such as downregulation of ISL-1and HCN4, and ectopic-expression of Cx40 and Nkx2.5 at the SAN site [18]. Finally, overexpression of Shox2 in embryonic stem cells can induce the formation of pace-like cells and increase the frequency of physiological pacing [19].
The 4q25 site mutations on human chromosome 4 are associated familial atrial fibrillation and are independent risk factors for early-onset atrial fibrillation [20,21]. Patients carrying this mutation are more likely to have Fig. 2 The expression of genes and proteins related to cardiac pacemaker cells in four groups. a Schematic of the experimental protocol for RA in four groups. b-i The mRNA level of the pacemaker markers. j-o The protein expressions of the pacemaker markers. Data points present as means ± SEM (n = 1). *p < 0.05 compared to CHIR = 0 group. # p < 0.05 compared to CHIR = 3 group. & p < 0.05 compared to RA + CHIR = 0 group a psychogenic stroke and recurrence after atrial fibrillation ablation [22,23]. The Pitx2 gene is close to this site and is associated with atrial fibrillation [24,25]. Pitx2 is widely expressed in the left side of the heart, such as the left atrium, pulmonary vein, and left vena cava, which induces ectopic electrical activity [26]. Wang's study showed that Pitx2 affects two microRNAs (miR-17-92 and miR-106b-25) and indirectly regulates the expression of Shox2 and TBX3, thus causing sinus node dysfunction and increasing the susceptibility to atrial fibrillation [27].
TBX5 is involved in the formation of the entire heart (cardiac conduction system, atrium, ventricles, etc.).
Initially, TBX5 is widely expressed in the posterior side of the first and the second heart fields of the heart at the embryonic stage [28]. TBX5 protein can act on the upstream regulatory sites of many transcription factors, such as TBX3, Shox2 and Pitx2 [29]. Moreover, Zhang et al. found that the Hedgehog-TBX5-HCN4 pathway regulates embryonic development of the atrial nodes of the heart, and Raghunathan's study showed that SHT5 (the combined overexpression of three factors, including Shox2, HCN2, and TBX5) can increase the proportion of pacing cells [30]. Additionally, TBX5 + Nkx2.5 − cells can label cardiac sinus node cells, while TBX5 + Nkx2.5 + cells The pacemaker cells rate marked by both Shox2 positive and Nkx2.5 negative (Shox2 + Nkx2.5. − ) in four groups detected by flow cytometry (n = 4). c Cells differentiated from four groups were stained for cTNT (green), Shox2 (red) and DAPI (blue). Magnification, × 200, n = 4 mark the first heart field [31]. Finally, TBX5 is found in the cardiac conduction system in adulthood [32]. TBX3 and TBX5 can jointly mark the cardiac conduction system, and TBX5 knockout mice show abnormal heart rates, such as sinus arrest, AV block and ventricular tachycardia [31,33]. Figure 8 shows possibly potential biological mechanisms of our article: on the one hand, RA can augment the expression of TBX5, Shox2 and HCN4. As mentioned above, TBX5 acts on the upstream regulatory sites of Shox2. On the other hand, RA can also reduce Pitx2 and miR-106b, two upstream inhibitors of Shox2 during the development of SAN. Thus, we inferred that RA-TBX5/Pitx2-Shox2-HCN4 pathway is a possible mechanism by which RA signaling initiates pacemaker differentiation. In addition, we found that RA had a negligible effect on the mRNA level and protein expression of ISL-1 and TBX3, which are mainly enhanced by Wnt signaling. Therefore, we combined both Wnt signaling and RA signaling to generate sinus node-like cells because of their synergistic effects. Moreover, previous studies have shown that BMP signaling specifies cardiac mesoderm toward the pacemaker cells fate, and our study disclosed RA can increase the mRNA level of BMP4 significantly. Finally, RA and CHIR together can slightly regulate TBX18 and Nkx2.5.
The indispensable role of the Wnt signaling pathway in SAN development was shown previously. Using tracer technology, researchers found that the Wnt signaling pathway is continuously expressed in the "third heart field," which is conducive to stage 4 automatic depolarization [34]. Researchers injected Wnt8c into the first and second heart regions and found that the speed of automatic depolarization and pulsating heart rate in stage 4 increased significantly. In addition, the classical wnt5b can directly induce Nkx2.5 + cells to sinus node-like cells [7]. Finally, another study suggested that the classical Wnt signaling pathway promotes the differentiation of cardiac precursor cells into sinus node cells, while the non-classical pathway regulates the formation of working cardiomyocytes [35,36].
Due to the limitations of current technology and the lack of specific markers for SAN, the source of SAN is not clear, although the current view is that the sinus node cells are derived from the right venous sinus [37][38][39]. Structurally and functionally, SAN is a heterogeneous tissue, i.e., it is a synthesis of many different subpopulations of cells. The sinus node consists of a large "head" within the right jugular vein myocardium and a "tail" distributed along the terminal crest [40]. TBX18 plays a decisive role in the "head" structure, which has a small effect on the electrophysiological function of the sinus node. In contrast, a Tbx3-deficient or Shox2-deficient embryo can Fig. 4 The frequency of beats and the expression of ionic channels in four groups. a The frequency of beats in four groups. b-h The mRNA level of ionic channels in four groups. Data points present as means ± SEM (n = 5). *p < 0.05 compared to CHIR = 0 group. # p < 0.05 compared to CHIR = 3 groups. & p < 0.05 compared to RA + CHIR = 0 group form a morphologically normal SAN but a slowed heart rate.
Based on the above results, researchers have used several small molecules to induce sinus node-like cells. Considering that the differentiation process is complex, it is necessary to seek a simpler and more effective differentiation scheme. Additionally, whether sinus node cells induced by those processes belong to the same subtype needs to be elucidated. Ren et al. found that Wnt signaling may initiate pacemaker differentiation by interacting with both ISL-1 and Tbx18 enhancers [7]. Protze et al. showed that BMP can increase the proportion of SANLPCs, while RA enhances the pacemaker phenotype of SANLPCs [8]. Moreover, the "SAN-linked cells" expressing Shox2 + Nkx2-5 + in the transition zone play a crucial role in sinus node cells and the excitation of atrial cells excitation [41]. The results of our flow cytometry assay showed that the proportion of Shox2 + Nkx2-5 + cells increased in RA + CHIR = 0 group, while in RA + CHIR = 3 group, Shox2 + Nkx2-5 − cells increased considerably. Wnt signaling pathway promotes the formation of structural TBX18 + "head" sinus node cells, RA increases the functional "tail" cells, and the combination of both might better mimic the formation of SAN. This further explains the decrease in the beat frequency induced by the overstimulation of Wnt signaling pathway, while the addition of RA can increase and restore the frequency of beats.
RA plays a pleiotropic role during heart development at different concentrations [42]. RA mainly induces . e Current density-voltage relationships in four groups (n = 28). f AP was detected in four groups. g APA, MDP level, 0-phase depolarization rate (dV/dt) and APD of four groups.*p < 0.05 compared to CHIR = 0 groups. APA, action potential amplitude; action potential duration; MDP, maximum diastolic potential epicardial cells at 1-4 µmol/L [43,44], 0.5-1 µmol/L of RA intervenes the mesoderm to increase the proportion of atrial myocytes [45]. High concentrations of RA (1-100 µmol/L) induce the production and maturation of ventricular muscle cells [46,47]. Liu and Hou conducted certain experiments [5,6] and found that RA inhibitors (BMS189453) can increase the proportion of sinus node cells, which appears to contradict our results at first glance. However, a closer analysis shows that BMS453 is an agonist of RARβ receptors, and an inhibitor of RARα as well as RARγ receptors. RA is pleiotropic as a broad agonist of RAR receptors and RXR receptors. After combining our results with those of Liu and Hou, we preliminarily inferred that the promotion of sinus node cell production facilitated by RA might partly be due to the effect on RARβ receptors. Some studies have shown that RA acts upstream of the transcription factor TBX5 and directly regulates the expression of TBX5 [48]. Exogenous RA-induced embryonic stem cells can increase transcriptional levels of Shox2 and ISL-1. These experiments showed that RA plays a key role in the development of SAN, and based on the results of Protze et al., Liu et al. as well as Hou et al., we focused our attention on investigating the RA signaling pathway. In the future, we aim to compare the different effects between BMS189453 and RA or explore the synergistic interaction with other signaling pathways.
The most prominent difference between our experiment and Protze et al. was the time of the intervention. Our experiment initially followed the RA intervention time suggested by the study of Protze et al., but we obtained atrial myocytes following that protocol. Their intervention point was mesoderm cells, but our intervention time was heart mesoderm cells and heart precursor cells, which was consistent with the study by Liu and Hou. Actually, the time of RA intervention in our experiments was consistent with the expression of RA in SAN embryonic development. In the E7.5-E8.0 stage, Raldh2 is initially expressed in the lateral side plate mesoderm, and then came to the pSHF cardiac precursor cells, limited to the sinoatrial node and the atrial area [49]. In order to exclude the specificity of the U1 cells, we induced another cell line (H9) following our protocol, and the results showed that RA can also increase the expression of Shox2 and reduce the mRNA level of Nkx2.5 (Additional file 7: Figure S3). Compared to the proportion of SANLPCs (35%) obtained by Protze et al., our proportion of SANLPCs (21.3%) was lower, but our protocol was more concise. Furthermore, we combined the use of RA in our study with other factors Fig. 6 The biological mechanism by which RA promotes pacing cell differentiation. a-e The mRNA level of the pacemaker markers on D20. f The mRNA level of TBX5 on D0-D19. g-h The protein expression of TBX5. Data points present as means ± SEM (n = 5). *p < 0.05 compared to CHIR = 0 group. # p < 0.05 compared to CHIR = 3 group. & p < 0.005 compared to RA + CHIR = 0 group (such as FGFi, BMP4) to determine whether the proportion of SANLPCs can be further increased.
Our experiments had three main limitations: Firstly, the administration of excess RA during pregnancy can have a teratogenic effect on the embryo. RA has a twosided effect on human gene regulation. Thus, whether it has some potentially detrimental effects (such as carcinogenicity) remains undetermined [50,51]. Secondly, previous studies have demonstrated that Shox2 acts as an upstream transcription factor of ISL-1 as well as TBX3. According to our results, RA improved the expression of TBX5 and Shox2 but did not affect TBX3 and ISL-1, the reason is not known. That is to say, the biological mechanism by which RA induces sinus node-like cells requires further investigation. Finally, due to the limitations of experimental conditions and taking into account the immunosuppressive effect between species, our experiment ended at the cell stage, and did not include the effects of the administration of RA on animals. Therefore, several animal-related problems should be considered in the follow-up experiments, for example, how to avoid xeno-immune rejection, whether RA has an operative effect in vivo and so on.

Conclusions
In conclusion, in this study, we found that RA with Wnt signaling can induce hiPSCs to differentiate into sinus node-like cells, increase the proportion of CTNT + Shox2 + Nkx2.5 − cells, and rise the frequency of beats. The underlying biological mechanism by which such changes occur might include cooperation with the Wnt signaling pathway to regulate pacing-related transcription factors, such as Shox2, TBX5, Nkx2.5, TBX18 and so on.